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Abstract 

We use the finite lattice method to count the number of punctured staircase and 
self-avoiding polygons with up to three holes on the square lattice. New or radically 
extended series have been derived for both the perimeter and area generating func- 
tions. We show that the critical point is unchanged by a finite number of punctures, 
and that the critical exponent increases by a fixed amount for each puncture. The 
increase is 1.5 per puncture when enumerating by perimeter and 1.0 when enumer- 
ating by area. A refined estimate of the connective constant for polygons by area is 
given. A similar set of results is obtained for finitely punctured polyominoes. The 
exponent increase is proved to be 1.0 per puncture for polyominoes. 

1 Introduction 

A self- avoiding polygon (SAP) can be denned as a walk on a lattice which returns to 
the origin and has no other self-intersections. Alternatively we can define a SAP as 
a connected sub-graph (of a lattice) whose vertices are of degree 2. The history and 



significance of this problem is nicely discussed in ||15| . Generally SAP are considered 
distinct up to a translation, so if there are p n SAP of length n there are 2np n returns to 
the origin (the factor of two arising from the two possible directions in which the loop 
can be traveled). Staircase polygons are a special case of SAP. A staircase polygon can 
be viewed as the intersection of two directed walks starting at the origin, taking steps 
either to the right or down and terminating once the walks meet. We define a punctured 
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SAP as a SAP with one or more holes, with the perimeter of each hole itself being a 
SAP. In other words a punctured SAP is a SAP enclosing one or more SAP. To avoid any 
possible confusion in our definition of punctured polygons, note that the punctures are 
disjoint — there are no degree four vertices in the punctured polygons we are considering. 
A similar definition can be given for punctured staircase polygons. We have shown an 
example of each case in figure [I]. The two principal questions one can ask are "how many 
polygons, distinct up to a translation, are there of perimeter 2n (including the perimeter 
of the holes) with k punctures?" and "how many polygons, distinct up to a translation, 
are there of area n with k punctures?" 

A polyomino is defined as the union of connected cells, where a cell is a unit square 
with 4 edges (and its interior). Two cells are said to be joined if they share a common 
edge, and are said to be connected if there exists a sequence of cells joining the two cells 
such that successive pairs of cells are joined. A fc-punctured polyomino is a polyomino 
with k holes. Unlike the situation for punctured polygons, degree 4 vertices are permitted, 
but two holes meeting only at a vertex do count as two holes, rather than one, as they 
are not joined. Punctured polygons are a proper subset of punctured polyominoes. The 
difference between them is that punctured polygons do not include those polyominoes 
where there are 2 cells touching each other only at a vertex. An example is shown in 
figure ||. 

Reverting to A;-punctured polygons, for k = we have the simpler — but still unsolved 
- problem of the number of polygons of perimeter In or area n. Both these problems 
have been extensively studied for several decades. The most recent result for polygon 
perimeters appears to be |17j where polygons of perimeter up to 90 steps are given. In 
that paper our analysis of the polygon perimeter generating function led us to conclude 
that 

p(°>(x) = ]>>Sx" ~ A®{x) + B^(x)(l - fi 2 x) 2 - a , 

n 

where p^l * s the number of SAP of perimeter 2n, distinct up to translations. The analysis 
in [17] yielded a very accurate estimate for the connective constant \x = 2.63815853034(10) 
and confirmed the theoretical prediction a = 1/2 [|19|. Furthermore, we obtained es- 
timates for the critical amplitudes A^°'(x c ) ~ 0.036 and B^°\x c ) ~ 0.234913, where 
x c = l/fi 2 . We also concluded that there was no evidence for a non-analytic correction- 
to-scaling exponent, so that we expect the asymptotic form of the coefficients to behave 
as: 

P2I ~ H 2n n~^[ai + a 2 /n + a 3 /n 2 + a 4 /n 3 H ]. 

The connective constant fi is of course the same as that for self-avoiding walks on the 
same lattice P~5"|. 



For polygon areas the most recent published work appears to be || in which the first 
20 terms of the area generating function were given and analysed. In that work it was 
found that 

A®(y) = J>£V ~ C {0) (y) + D {0 \y) log(l - «y), 
n 

where is the number of SAP of area n, k fa 3.97087, and the amplitudes and 
were not estimated. Recently this series has been extended to 26 terms |20| , but in 



the present work we have devised a new and exponentially faster algorithm, as a result of 
which we have extended the series to 42 terms, and we present an analysis of this longer 
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series. The connective constant k is found to be slightly smaller than that for the related 
problem of polyominoes [|T^ . 

In the following, we refer to the boundary of a polygon and its interior as a disc and 
so we will be discussing punctured discs. An unpunctured disc is a SAP. 

For punctured discs, the basic problem is, analogously, the calculation of the generating 
functions 



P( fe )( x ) = J2P2n xU ~ B{k \ X ) + C(k) faX 1 " (/i (fc) )^) 2 ~ Qfe , (1) 
n 

and 

A«(y) = £ «£V ~ ^ (fc) (y) + £ (fc) (2/)(l - /c«y)-*, (2) 

n 

where the superscript k refers to the number of holes, or punctures. From the generating 
function, one wishes to deduce the asymptotic behaviour, believed to be as shown on the 
r.h.s. of the above equations. The major problem to be investigated is how the behaviour 
of (x) and (y) changes as k is increased. Previous work |25|, ^] has been confined 
to the study of punctured SAP by area. There it was proved that K^ k ' = k^ ' = k, and 
that if the exponent exists, Pk = Po + k. These results apply more generally to punctured 
surfaces, but in this work we are confining ourselves to two dimensions. As far as we 
are aware, there has been no previous work on the problem of the perimeter generating 
function of punctured discs. 

The problem is interesting for several reasons. The effect of a change in geometry 
is a much studied topic in lattice statistics, and our study of the change in perimeter 
exponent with punctures seems to be entirely new. It has only been possible by the 
algorithms we have designed and implemented, which are exponentially faster than pre- 
existing algorithms. A number of related models have been studied previously, such as 



c- animals [181 M and the behaviour of prime knots in polygons [22]. c-animals are lattice 



animals with exactly c cycles. In [M] it was proved that if the number of such animals 



a n (c) ~ A"n c as n — > oo, c fixed, then 9 C = 9q + c provided that 9q exists. It had been 



previously proved |28| that A c = Ao- The change in connective constant of c-animals as 
the number of cycles per vertex changes from zero to non-zero is discussed — among 
other results — in [18|. Similarly, in a numerical study of knotted polymers p2 , it was 



conjectured that the exponent a depends on the number of prime knots n p that arise in 
the knot decomposition of a given SAP via the relation a(n p ) — a(0) + n p . 

Furthermore, there is considerable pedagogical connection between some of these pre- 
vious studies and our work here, and also between our work and the study of branched 
polymers. To sketch this connection, we first remark that the number of polyominoes, 
also called lattice animals, is just the number of (strongly embedded) site animals. This 
connection is readily seen by placing a site at the centre of every cell of the polyomino, 
and joining those sites corresponding to joined cells by a bond. In this way, every lattice 
animal is mapped to a distinct site animal and vice versa. All other models we are consid- 
ering map similarly to a different subset of strongly embedded site animals. For example, 
SAP map to site animals whose only cycles are 4-cycles (which may be isolated or joined). 
Punctured polygons map to site animals with larger cycles, as well as 4-cycles, whereas 
punctured polyominoes correspond to site animals with more complex restrictions. Thus 
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this study also complements the earlier studies of c-animals. In those studies, the vari- 
ation of exponent with the number of cycles is considered, whereas in this study we are 
varying the types of allowed cycles. 

In order to study these and related systems, when an exact solution can't be found 
one has to resort to numerical methods. For many problems the method of series expan- 
sions is by far the most powerful method of approximation. For other problems Monte 
Carlo methods are superior. For the analysis of P ( - k \x) and A^ k \y), series analysis is 
undoubtedly the most appropriate choice. This method consists of calculating the first 
few coefficients in the expansion of the generating function. Given such a series, using the 
numerical technique known as differential approximants ||13| , highly accurate estimates 
can frequently be obtained for the critical point and exponents, as well as the location 
and critical exponents of possible non-physical singularities. Other numerical methods 
are discussed in |T7|], and those used in the present study are described more fully below. 



In the next section we will describe the finite lattice method for enumerating punctured 
polygons. In Section 3 we prove the invariance of fx^ h ' as k changes, and give an heuristic 
argument for the exponent shift with k. The results of the analysis of the series are 
presented in Section £|. Results analogous to those known for punctured polygons by area 
are proved for punctured polyominoes in the Appendix. 



2 Enumeration of punctured polygons 

2.1 Enumeration of punctured self-avoiding polygons 

The method used to enumerate punctured self-avoiding polygons on the square lattice 
is a generalisation of the method devised by Enting |7j for the enumeration of ordinary 
SAP. In the following we first describe the original method in some detail and show how 
simple it is to generalize the method to the enumeration of punctured polygons. The first 
terms in the series for the generating function can be calculated using transfer matrix 
techniques to count the number of polygons in rectangles W + 1 edges wide and L + 1 
edges long. The transfer matrix technique involves drawing a line through the rectangle 
intersecting a set of W + 2 edges. For each configuration of occupied or empty edges along 
the intersection we maintain a (perimeter) generating function for loops to the left of the 
line cutting the intersection in that particular pattern. Polygons in a given rectangle are 
enumerated by moving the intersection so as to add one vertex at a time, as shown in 
figure ^. Since the loops are non-intersecting, each configuration can be represented by 
an ordered set of edge states {rii}, where 



{0 empty edge, 
1 lower part of loop closed to the left, 
2 upper part of loop closed to the left. 

Configurations are read from the bottom to the top. So the configuration along the 
intersection of the polygon in figure ^| is {0112122}. In passing it is worth noting that 
there are some major restrictions on the possible configurations. Firstly, since all loop- 
ends are connected to the left of the intersection, every lower loop end must have a 
corresponding upper end, and it is therefore clear that the total number of 'l's is equal 
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to the total number of '2's. Secondly, as we look through the configuration starting from 
the bottom the number of 'l's is never smaller than the number of '2's. 

In Table [3] we have listed the possible local 'input' states and the 'output' states which 
arise as the kink in the intersection is propagated by one step. Some of these update 
rules are illustrated further in Fig|4j. The first panel represents the input states '10' and 
'01' and the possible output states are also '01' and '10'. The second panel represents 
the input state '11' as part of the configuration {01122}. In this case we connect the two 
loop ends, but in doing so we see that the upper part of the second loop before the move 
becomes the lower part of the one remaining loop after the move. That is the configuration 
{01122} becomes {00012}. This relabelling of the other loop-end when connecting two 
'l's (or two '2's) is denoted by over-lining in Table [1|. In general there could be more loops 
nested in between the two 'l's and the corresponding '2' at the other end of the loop. Say 
for instance we had the configuration {11121222} and connected the first two 'l's then 
the new configuration of unconnected loop ends would be {00121212} (drawing a little 
figure makes this quite clear). The general rule for the relabelling is as follows: When 
connecting two 'l's ('2's) we work upward (downward) in the configuration, counting the 
number of 'l's and '2's we pass until the number of '2's ('l's) exceeds the number of 'l's 
('2's). This '2' ('1') is the other end of the inner loop and it should now be changed 
to a '1' ('2'), thus becoming the lower (upper) end of the outer loop (again drawing a 
few pictures should make this clearer). The weights corresponding to these configuration 
transformations are simply calculated by counting the number of steps which have been 
added to the polygon. Note that the input state '12' is special because connecting the two 
ends results in a closed loop, so this is only allowed if there are no other loops cut by the 
intersection and the result is a valid polygon, which is then accumulated in the total count 
for that particular length. Failure to observe this restriction would result in graphs with 
disconnected components, either one polygon over another or a polygon within another 
(this latter case is of course of interest when we wish to enumerate punctured polygons). 
This is illustrated in figure |5] where we show the possible ways a pair of loops can be 
placed relative to one another (first panels), how the loops can be connected to produce a 
valid configuration (middle panels) and the ways of connecting loops that lead to invalid 
graphs containing disconnected components (last panels). In this figure the invalid SAP 
in the last panel on the bottom could result in a valid punctured SAP. Note also that 
from the input state '00' we can produce the output state '12' only if there are other loops 
crossing the intersection (otherwise we would produce disconnected polygons sitting side 
by side). We refer the interested reader to j7| || for further details regarding the encoding 
and relabelling of configurations. 

Due to the obvious symmetry of the lattice one need only consider rectangles with 
L > W. In the original approach valid polygons were required to span the enclosing 
rectangle in the lengthwise direction. So it is clear that polygons with projection on the 
2/-axis < W, that is polygons which are narrower than the width of the rectangle, are 
counted many times. It is however easy to obtain the polygons of width exactly W and 
length exactly L from this enumeration M. Any polygon spanning such a rectangle has 
a perimeter of length at least 2(W + L). By adding the contributions from all rectangles 
of width W < W max (where the choice of W max depends on available computational 
resources) and length W < L < 2W max — W + 1, with contributions from rectangles with 
L > W counted twice, the number of polygons per vertex of an infinite lattice is obtained 
correctly up to perimeter 4W max + 2. 
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With the original algorithm the number of configurations required as W mSLX increased 
grew asymptotically as 3 Wn1 ^ |12|]. In a recent improvement of the algorithm QIBL |T7 



valid polygons were required to span the rectangle in both directions. In other words we 
directly enumerate polygons of width exactly W and length L. For each configuration 
of partially completed polygons we keep track of the current minimum number of steps 
N cm that have been inserted to the left of the intersection and we calculate the minimum 
number of additional steps iV a dd required to produce a valid polygon that spans a rectangle 
of size at least W x W. If the sum N cur + iV add > 4W niax + 2 the partial generating 
function for that configuration was discarded because it would make no contribution to 
the polygon count up to the perimeter lengths we were trying to obtain. Numerical 
evidence indicated that the computational complexity was reduced significantly. While 
the number of configurations still grew exponentially as A H/max the value of A was reduced 
from A = 3 to A ~ 2 with the improved algorithm. Furthermore, for any W we know that 
contributions will start at 4W since the smallest polygons have to span aWxW rectangle, 
so for each configuration we need only retain 4(W max — W) + 2 terms of the generating 
functions while in the original algorithm contributions started at 2W + 2 because the 
polygons were required to span only the length-wise direction. 

The generalization to enumeration of punctured polygons is obtained by noting that 
a closed loop is formed whenever we connect a 1-edge to a 2-edge immediately above. If 
these two edges were the only loop ends in the intersection we would have formed a valid 
polygon. In other cases we need to ensure that the resulting polygon is a valid punctured 
polygon. That is, we must ensure that the separate SAP just formed will be completely 
enclosed within a larger polygon. So we wish to avoid forming separable polygons, as 
shown in the upper right panel of figure [| and 'holes within holes'. As it turns out the 
rule for a valid '12' closure is simply that we can connect the two loop ends provided 
there is an odd number of loop ends below the loop being closed. To see this consider 
that as we go through a configuration we note that each time we pass a loop end we go 
from the outside of the polygon to the inside and visa versa. So in a configuration all 
lattice cells between the first and second loop ends will lie inside the finished polygon, 
lattice cells between the second and third loop end will lie outside the polygon, and so 
on. Thus we see that by closing a loop which has an odd number of loop ends below it we 
are closing off a part of the lattice which will lie inside the finished polygon. In particular 
we see that we avoid the situation shown in the upper right panel of figure |5] with graphs 
containing disconnected pieces one over another. Likewise we avoid creating graphs with 
disconnected pieces sitting side by side. 

We also avoid forming holes within holes because closing a loop around the hole would 
be prohibited since there would be an even number of loop ends below the 'l'-edge (except 
of course when forming a completed punctured polygon). Let us look at the possible 
edge-configurations around a puncture in some detail. First look at the configuration 
{. . . 1122 . . . } (where the . . . are any edge configurations with an even number of edges 
making the total configuration of edges a valid intersection). The outer '12'-edges can't 
be connected (the number of edges below the '1' is even) so we have to connect two other 
edges on either side of the hole diminishing the number of edges by two and possibly 
changing the edge labels on either side of the hole (in which case we end up with one of the 
subsequent cases). In the top panel of figure || we show how the simplest interesting case 
{11112222} leads to only valid punctured polygons. Secondly, in the case {. . . 1121 . . . }, 
we can connect the two '1 '-edges forming a partial loop enclosing the hole. In doing so 
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the matching '2'-edge of the second 'l'-edge is changed to a 'l'-edge because it now is the 
new lower edge of the larger loop formed by connecting the two original 'l'-edges. The 
important thing to notice is that the number of edges below this new 'l'-edge is even 
so we cannot connect it to a '2 '-edge immediately above and we therefore do not form 
a hole within a hole. In the middle panel of figure ^ we show the simplest interesting 
case {1111212222} and demonstrate that it leads to valid punctured polygons. Thirdly, 
the configuration {...2121...} is not very interesting since connecting the '21'-edges 
in front of the hole does not result in a loop partially enclosing the hole. Finally, the 
configuration {. . . 2122 . . . } is obviously merely the mirror image of the second case. Note 
that this covers all cases since more complicated configurations merely would correspond 
to more convoluted loop structures. 

In this work we use the generalisation of the algorithm of |17j to count the number of 
punctured SAPs with up to 3 holes. Obviously, the smallest hole we can make in a SAP 
has perimeter 4 so the number of punctured SAPs is obtained correctly up to perimeter 
4W max + 2 + Ak. 

The algorithm used for the enumeration of punctured SAPs by area is a simple varia- 
tion of the algorithm described above. The encoding of configurations along the intersec- 
tion and the transformation of these configurations as the intersection is moved remain 
the same. The only change is that the weights are different. In order to count the enclosed 
area we proceed as follows: A unit of area may be added as the kink is moved to a new 
lattice cell (in figure |3| this is the cell in which the dotted lines meet). Whether or not a 
unit of area is added is determined by whether or not this lattice cell is inside or outside 
the polygon. But we already know from the arguments given above that a lattice cell is 
inside the polygon if the number of loop ends below the cell is odd (note that any loop end 
along the vertical edge cutting the horizontal part of the kink is included in the count). 
So in the case shown in figure |^ the lattice cell to which the kink is moved lies outside 
the polygon (there are two loop ends below the kink) and no unit of area is added. Note 
that a unit of area may be added for any given output configuration. In this case the area 
generating function is obtained correctly to 2W max + 3k. The factor 3k arise since it takes 
3 lattice cells to completely surround the simplest puncture, which is just a single cell. 

As far as we are aware, series for punctured polygons have not been derived previously. 
In such circumstances it is even more important than usual to undertake careful tests to 
ensure that the series are correct. To this end, a second algorithm was implemented, 
by one of the authors, to independently evaluate the series coefficients. The complete 
agreement we obtained between the two data sets reassures us as to the correctness of 
our results. 



2.2 Enumeration of punctured staircase polygons 

The enumeration of punctured staircase polygons is much simpler. In fact as we shall 
demonstrate it is a problem for which the computational complexity grows only as a 
polynomial in the number of terms. As stated in the Introduction we can think of stair- 
case polygons as consisting of two non-intersecting directed walks on the square lattice. 
Punctured staircase polygons naturally will have more than two walks, in fact up to 2k + 2 
walkers can be present in any given column. Due to the restrictions on staircase polygons 
it follows that a punctured staircase polygon is formed by requiring that the walks be 
directed and that any two walkers starting at one point join each other later without 
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intersecting other walks. This can be encoded in a transfer matrix calculation as follows: 
Again we count the number of punctured polygons in rectangles as before and draw a line 
through the lattice intersecting the W + 2 edges. In this case we need only specify whether 
or not each edge is part of the polygon or not, so each configuration can be represented 
by an ordered set of edge states {rii}, where 

J empty edge, 

[1 part of loop closed to the left. 

This uniquely specifies the configuration because the first occupied edge is connected to 
the last occupied edge and any edges in between are paired, e.g., the second and third 
(fourth and fifth and so on) edges form a loop to the left and has to be connected to each 
other later on. The rules for updating the configurations are as follows: From the local 
input state '00' we can always get the output '00', and the output '11' provided there is 
an odd number of edges below the kink and also that the edge directly below the kink is 
empty. The rules for the output '11' ensure that the new pair of walkers lie within the 
enclosing staircase polygon but not inside an internal staircase polygon (thus preventing 
holes within holes), and that the lower walk does not intersect other walks. From the 
local input state '01' and '10' we can always produce the output '01', and the output '10' 
provided either that the edge directly below the kink is empty or there is an even number 
of edges below the kink. These rules ensure that the walkers do not intersect other walks 
except when we close a valid loop. Finally from the local input state '11' we always get 
the output state '00'. The weights associated with these updates are obtained in the same 
way as for the SAP enumeration, whether the enumeration is done by perimeter or area. 

As in the previous case we calculate the number of punctured staircase polygons 
spanning the rectangles. Adding the contributions from all rectangles of width W < W mSuX 
and length W < L < 2W ma _ x — W + 1 the number of punctured staircase polygons is 
obtained correctly up to perimeter 4W / max + 2 + 4k. Note that for fixed k the maximum 
number of configurations Nq grows as a polynomial in W max 

3=0 V 7 

So in this case the algorithm is of polynomial complexity. 



3 Expected behaviour 

As mentioned in the introduction, the problem of enumerating SAP by area has been 
extensively studied || ED|. It has been shown that = k^°\ and that if the exponents 
exist, Pk = Po + k ||26|, ^7|. That is to say, for k finite, the connective constant for k- 
punctured discs, by area, is the same as that for unpunctured discs, while the critical 
exponent increases by 1 for each puncture. We repeat that the punctures are disjoint. 

As far as we are aware, there has been no previous work on the problem of the perimeter 
generating function of punctured SAP. We first give a proof that a 1-punctured SAP on 
the square lattice has the same connective constant as an unpunctured SAP, and indicate 
how that proof can be generalised to £;-punctured SAP. 
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Before stating and proving the relevant theorem, we require certain preliminary results. 
For an (unpunctured) SAP of perimeter 2n one knows, |I5| eqn. (7.101), that 



exp[-bVn]v 2n /4n<p ( £<fi 2n , (3) 

where b is a constant. Further, it is obvious that a polygon of perimeter 2n has maximum 
area n 2 /4, which occurs when the shape is an n/2 x n/2 square. Thus the area of a 
polygon of perimeter 2n, denoted Aim satisfies Am < n 2 /A. 

Now consider 1-punctured polygons of total perimeter 2n, with inner perimeter 2m 
and outer perimeter 2(n — m). It is possible for the inner polygon to be Hamiltonian, in 
which case its perimeter is greater than that of the surrounding polygon. Indeed, a square 
polygon of side 2n + 1, and hence of perimeter 8n + 4 can contain an internal polygon of 
perimeter as large as An 2 . Hence the semi-perimeter of the inner polygon, m, can range 
from a minimum value of 2 to a maximum value of n — 2y2n + 2 < n — y/n for n > 1. 

With these preliminaries, we can now state and prove the following theorem: 
Theorem: lim^oo j- log p^ = log/i, where fi is the same constant as appears in the 
corresponding limit for unpunctured SAP. 

Proof: 1-punctured polygons of total perimeter 2n are constructed by placing polygons P 
of perimeter 2m inside polygons Q of perimeter 2n — 2m. Let w(P, Q) denote the number 
of ways of placing polygon P inside polygon Q, (which of course depends on both P and 
Q.) Then 

n—y/n 

m=2 P Q 

We bound this summand by the product of three factors. The first two factors are the 
number of polygons of perimeter 2m and 2n — 2m respectively. The third is the number 
of ways the smaller polygon can be placed inside the surrounding polygon and is clearly 
less than or equal to the area of the surrounding polygon. Explicitly, 

pS < ( n " ™)W (n ~ m) /4 < yu 2 74 ^(n - mf < n 3 fi 2n . 

m=2 m=2 

To obtain a lower bound, consider a 3 x 3 square polygon with a unit square hole 
at its centre. This unique realisation of can be uniquely concatenated with each 
unpunctured polygon by joining them at a solitary specified edge, and then deleting that 
edge. For each unpunctured polygon we take the set of left-most vertical edges and choose 
the bottom edge from this set, and we choose the right-most, top-most vertical edge of p^ 
(a similar operation is shown for polyominoes in figure AT). The concatenation operation 



then gives, for each of the P2n-i4 unpunctured polygons, a unique member of the set of 
1-punctured polygons in which the puncture is a single cell. Thus we have p^-u — Pzn- 
From the above equations we thus obtain 

e W [-bVn~7]^ 2n - u /A(n - 7) < p^-u < pH < n^ 2n 

for n > 8. The theorem then follows immediately on taking logarithms, dividing through 
by 2n and taking the limit as n — ► oo. 
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This proof can clearly be extended to two-punctured polygons, then to three-punctured 
polygons etc., by concatenating unpunctured polygons with minimal two-punctured, three- 
punctured etc. polygons. In the appendix we prove that fc-punctured polyominoes have 
the same growth constant as unpunctured SAP by area. 

We have been unable to prove a result analogous to the result for the exponent of 
punctured SAP by area, but give an argument that depends on certain assumptions that 
are generally accepted, though not proved. In the case of staircase polygons however our 
assumptions have been proved, and so our result will be rigorously true. 

The key results we need are those obtained in || to the effect that the mean area of 
SAP of perimeter 2n is proportional to n 1 ' 5 . This is true both for SAP and for staircase 
polygons, and in the latter case it has been proved. More precisely, we need the following 
result. There exists constants Di and D 2 such that the mean area A 2n of polygons of 
perimeter 2n satisfies 

D 1 n$<A 2n = ^^-<D 2 n$, 

P2n 

where the sum is taken over all p 2n polygons Q of perimeter 2n. 

Further, there exists constants C\ and C 2 such that the number of polygons of perime- 
ter 2n satisfies 

CV"n"t <pg <C 2f i 2n n-l 

For staircase polygons, /i = 2 and the exponent is § instead of \ in the above equation. 

As above, let Q denote a polygon of perimeter 2n — 2m, with area Aq. The number 
of ways of placing a given polygon P of perimeter 2m inside Q is clearly less than Aq. 

Thus the number of ways of placing P inside any polygon of perimeter 2n — 2m is less 
than 

^2 Aq = p 2 l_ 2m A 2n _ 2m , (4) 
Q 

where the sum is over all polygons Q of perimeter 2n — 2m. 

Hence the number of ways of placing all of the pf^ polygons of perimeter 2m inside 
any polygon of perimeter 2n — 2m is less than p^p 2 \l_ 2m A 2n - 2m , and so 

p£ < E pSM-am^a™ < AAV"2- 5 S 2n , (5) 

m=2 

where 

S 2 n = E m ~V(™ - m )- 
m=2 

The last sum may be evaluated in a variety of ways. Using Maple we find that 

S 2n = 0.341../n + O(n _ 5). 

Thus we obtain the bound 

p2J < const. fi 2n /n. (6) 
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For staircase polygons the analogous calculation is slightly simpler, and we find 

S 2n ~l + ((3/2) + 0(n- l 2), 

where ((z) is the Riemann zeta function. 

To obtain a lower bound, we restrict the inner polygon to be of perimeter 4, that is, a 
unit square. The number of polygons punctured by a unit square clearly provides a lower 
bound to the number of 1-punctured polygons. A unit square can be placed anywhere in 
a polygon Q of perimeter 2n — 4 except on a boundary site. There are 2n — 4 boundary 
sites. The mean area of Q is greater than D\(n — 2)3, so we obtain the bound 

2n-4 

> P2n-A(D 1 (n - 2)3 - 2n + 4) > C x — M x (n - 2)3 - 2n) > const. ^ n jn. 

(2n-4)3 

(7) 

Combining the two bounds gives the result 

E lf i 2n /n<p$<E 2f i 2n /n, 

where E\ < E 2 are constants. Accepting the usual asymptotic form that is expected for 
such models, we conclude that 

p£ ~ const.fi 2n /n. (8) 

(For staircase polygons the analogous result is p^ ~ constA n .) Since the number of 
unpunctured polygons grows like /i n n _5//2 , we see that the exponent is predicted to increase 
by I as the result of a single puncture, while in the case of the area generating function, the 
exponent is found to increase only by 1. We show in the next section that this prediction 
is borne out by our numerical calculations. 

For punctured staircase polygons, a similar conclusion holds. That is, the connective 
constant is unaltered at \i sta ir = 2, but the exponent increases by 1.5 over its unpunctured 
counterpart when enumerating 1-punctured staircase polygons by perimeter. In fact, we 
have been able to calculate the generating function for staircase polygons with a single 
puncture of perimeter 4, and also with a single puncture of perimeter 6. The generating 
functions for these special cases are given below, and are precisely in accordance with the 
more general results given above. 

We also note that Cardy recently considered the problem of the number of punc- 
tured SAP with k concentric, mutually self-avoiding SAPs surrounding a fixed point of 
the dual lattice. Thus for k — 1 this corresponds to 1-punctured SAPs surrounding a 
fixed point. In that case Cardy finds || for the number of such configurations 6 2 n that 
^2n = /^ 2 " 2 n • The above calculation may be repeated for Cardy's problem. All that 
is required is to add a factor A m to the summand in eqn. ||, as the surrounded point 
can be anywhere inside the inner polygon. Thus we must multiply by the mean area of 
that polygon. The sum defining S 2n then becomes S 2n = X^m=2 m(J-m) ~ log n/n, in 
agreement with Cardy's result. (Our lower bound for this problem is too weak, but with 
more effort could be improved.) 
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4 Analysis of the series 



All the series we have investigated are characterised by coefficients that grow exponen- 
tially, with sub-dominant term given by a critical exponent. The generic generating 
function behaviour is G(z) = ^2 n g n z n ~ A(z){l — az ) an d hence the coefficients of 
the generating function g n = [z n ]G(z) ~ A(l/a)/T(^)a n n^~ l . Generally speaking the ex- 
istence of the growth constant a has been proved, but except for exactly solvable models, 
such as staircase polygons, the existence of the critical exponent £ has only been con- 
jectured, though its existence has never been doubted. The radius of convergence of the 
generating function is usually given by the critical point, which is at z — 1/a. 

We principally used two methods to analyse all the series studied in this paper. 
Firstly, to obtain the singularity structure of the generating function we used the nu- 
merical method of differential approximants [pl|. In particular, we used this method to 
estimate the growth constant a and the critical exponent £. We were invariably able to 
conjecture an exact value for £, which was always an integer or half-integer for all the 
problems we investigated. Imposing this conjectured exponent permitted a refinement of 
the estimate of the growth constant — providing so-called biased estimates. 

Once the exact value of the exponent was conjectured, and the growth constant ac- 
curately estimated, we turned our attention to the "fine structure" of the asymptotic 
form of the coefficients, by fitting the coefficients to the assumed form g n = [z n ]G(z) ~ 
a n n^~ 1 ^2 i>Q Ci/n^' l \ If there is no non-analytic correction term, then f(i) = i, while a 
square- root correction term means f(i) = i/2. For all the series studied, only these two 
situations were encountered. 

In all cases, our procedure is to assume a particular form for f(i), and see how it 
fits the data. With the very long series we now have at our disposal, it is usually easy 
to see if the wrong assumption has been made — the sequence of amplitude estimates 
Cj either diverges to infinity or converges to zero. Once the correct assumption is made, 
convergence is usually rapid and obvious. A detailed demonstration of the method can 
be found in §, |T7|. 

As an example of the sort of results we obtained, we show sequences of estimates of the 
coefficients of the perimeter generating function of unpunctured SAP in Table |[ In that 
case we conjectured that f{i) = i. Because of the large amount of tabular data generated 
by the method, we have not given this level of detail for the many series investigated 
here. We show only the results for two series. For the others, we just give our assessment 
of the apparent convergence of the sequences q, and the estimated value of the limits. 
As the equations involved are linear, the method is easy to implement, and interested 
readers can readily generate the relevant data themselves. Some subtleties nevertheless 
exist. For example, for punctured staircase polygons, the perimeter generating function 
has two singularities on the circle of convergence, and so both must be taken into account. 
We discuss this in more detail in the relevant section below. 

As for the first stage of the analysis, the method of differential approximants, we pro- 
ceeded as follows: Estimates of the critical point and critical exponent were obtained by 
averaging values obtained from first order [L/N; M] and second order [L/N; M; K] inho- 
mogeneous differential approximants. For each order L of the inhomogeneous polynomial 
we averaged over those approximants to the series which used at least the first 80% - 90% 
of the terms of the series, and used approximants such that the difference between N, M, 
and K didn't exceed 2. These are therefore "diagonal" approximants. Some approximants 
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were excluded from the averages because the estimates were obviously spurious. The error 
quoted for these estimates reflects the spread (basically one standard deviation) among 
the approximants. Note that these error bounds should not be viewed as a measure of the 
true error as they cannot include possible systematic sources of error. However systematic 
error can also be taken into account in favourable situations, as for example in the case of 



SAP enumerated by perimeter [ IT] . Again, in the interests of space, we present only our 
results, and not the intermediate detail from which our estimates were made. An example 
in full detail for one of the series investigated in this study can be found in JH]]. We turn 
now to the analysis of all the series. 



4.1 Staircase polygons 

For (unpunctured) staircase polygons, the multi-variable width, height and area generat- 
ing function is known 0. As usual, we denote 

n— 1 

(a) n = l[(l-aq l ). 

i=0 

Further, denoting the first two g-Bessel functions as: 

(-l) n x n q( n i 1 ) 

n>0 

and 



(q)n(yq)n 
n ^ (g)n-i(yg)n-i(i -yq n ) : 

the perimeter and area generating function is simply 



Ji(x,y,q) = ^2 



P{x,y,q) = y-r-, r, 

Ja\x,y,q) 

where x (y) is the variable conjugate to the horizontal (vertical) semi-perimeter, while q 
counts the area. No analogous results are known for punctured staircase polygons, though 
the calculation for SAP [26|] by area can be carried over mutatis mutandis to prove that 



the area generating function for a punctured staircase polygon with k holes has the same 
radius of convergence as the area generating function for unpunctured staircase polygons. 
Further, the critical exponent increases by 1 for each puncture, and unlike the case for 
SAP, we can not only prove the existence of a critical exponent for unpunctured staircase 
polygons, but we know its value. Hence we know the leading term in the asymptotic 
expansion of the generating function by area for /c-punctured staircase polygons. 

For the expected behaviour of the perimeter generating function of punctured staircase 
polygons, the arguments of the preceding section apply directly. The radius of conver- 
gence, and hence the connective constant remains unchanged, and the argument given in 
the preceding section suggests that the critical exponent should increase by 1.5 for each 
puncture. The results of our analysis, presented below, bear this out. 
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4.1.1 Area generating function 

For unpunctured staircase polygons, the area generating function is given by 

By inspection, this has poles at the zeros of J (l,l,g). The nearest zero is at |L] 
1/q = i] = 2.30913859330, and there is a simple pole at that point. The next zero is well 
separated (at 1/q = A = 1.4435..) and so the asymptotic form of the generating function 
is 

A{q)~D/(l- m ) + E/{l-\q) + ... , 

and hence 

aW = [q n }A(q)~ v n (c + O((\/ V )r). 

Our analysis bears this out, and we estimate c = 0.12881579. 

For 1-punctured discs, our analysis, based on more than 100 series coefficients, con- 
vincingly suggests the following asymptotic form: 

A (1 \q) ~ - mf + E {1) /(1 - m) 1 ' 5 + F {l) /(l - V q) + ..., 

and hence 

W = [q n ]A {1 \q) ~ v n nJ2ci/n i/2 . 

i>0 

The sequences of amplitude estimates, assuming this asymptotic form, are shown in Ta- 
ble [| The apparent convergence of the amplitude estimates is, as explained above, our 
source of evidence for this asymptotic form. 

For 2-punctured discs, a similar analysis, based on some 86 series coefficients, convinc- 
ingly suggests the following asymptotic form: 

A®(q) ~ D®/(1 - r,qf + E^/(l - V q) 2 - 5 + ..., 
a W = [q n ]AW(q) ~ 77V ^ ajn 1 ' 2 . 

j>0 

For 3-punctured discs our analysis was based on an 89 term series. We found that the 
above pattern persists, so that the generating function has the following asymptotic form: 

A^\q) ~ D (3) /(l - mf + E {3) /(1 - r/g) 3 - 5 + . . . , 

and hence 

ag) = [q n }A^(q) ~ ^n^ai/rj' 2 . 

i>0 

Estimates of the various amplitudes defined above are shown in Table |j. 
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4.1.2 Perimeter generating function 

For unpunctured staircase polygons, the perimeter generating function (ignoring the dis- 
tinction between height and width) is given by 



„, . 1 - 2x - VI - 4x 

P(x) = , 

\ j 2 

which is, apart from suppression of the first term, the generating function for Catalan 
numbers. Hence 

= [x»]P(x) = ^f^ 2 ) - r/nl^/n*. 

^ ' ?>o 

The values of q follow immediately from the exact solution. They are given in Table |j. 

For 1-punctured discs, our analysis, again based on more than 100 series coefficients, 
was more equivocal than that of the 1-punctured area generating function. The method 
of differential approximants clearly identified a singularity at the known critical point, 
x c = 1/4, but almost all approximants had a double root, implying a confluent singularity. 
The leading exponent was estimated to be —1, implying a pole in the generating function, 
(and hence immediately lending support to our conjectured change in the critical exponent 
of 3/2 as a result of puncturing), but we were unable, from this method, to identify the 
confluent exponent. Further, a second singularity was identified, of the form const. (1 + 
x/x c ) 6 ' 5 . At this point we wish to remark on the close similarity between a recently solved 
model of polygons, the three-choice polygon model and 1-punctured staircase polygons. 
1-punctured staircase polygons can be thought of as being constructed from two three- 
choice polygons, with common edges deleted. This geometric similarity is borne out by the 
fact that the two models have identical singularity distributions, with even the exponents 
being the same at all non-physical singularities. 

For 1-punctured discs then, the perimeter generating function is expected to be of the 
following asymptotic form: 

pW( x ) = J2p ( £x n ~ B {l \x) + C w (x)(l - Ax)' 1 + D w (x)(l + Ax) 6 ' 5 , 

n 

and hence 

pCD = [x n }P {l \x) ~ 4 n ^ Ci /™ /W + (-4) n rr 7 - 5 ^d;M 

i>0 i>0 

where B^(x), C^(x) and (x) are assumed regular in the disc \Ax\ < 1. 

Assuming f(i) = i, which implies only analytic correction-to-scaling terms, gave un- 
satisfactory results. Notably, we observed that the estimates of the amplitude C\ in the 
above asymptotic form were steadily increasing, suggesting that the assumed form did 
not properly account for the correction-to-scaling terms. With f(i) = i/2 the amplitude 
estimates were much more stable. This then implies that the generating function in fact 
behaves as 

P (1) (x) ~ P (1) (x) + C {1) (x)(l - 4a;)- 1 + P (1) (x)(l - Ax)^ + D {1 \x)(l + Axf 5 . (9) 

The amplitude estimates, assuming this asymptotic form, are shown in Table (|. The 
apparent convergence of the amplitude estimates is, as explained above, our source of 
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evidence for this asymptotic form. We did not tabulate the non-physical amplitudes, 
as they are of little interest to our investigation. However they need to be included to 
stabilise estimates of the physical amplitudes. We mention in passing that do ~ 0.14. 

We also generated series for staircase polygons with two and three punctures, with 
perimeter 150 and 134 steps respectively. Our differential approximant analysis lent sup- 
port to our expectation that the critical exponent increases by 1.5 per puncture. This 
turned out to be true for the non-physical singularity at X — ~7 clS well. A similar anal- 
ysis to that for the one-punctured polygons strongly supported the analogous asymptotic 
forms, 

p( fc ) ( x ) = P2nX n ~ B {k) (x) + C {k) (x) (1 - 4x) a5 - L5fc + D w (x) (1 + 4x) 8 - 1 - 5fe (10) 

n 

for k > 0. Hence the asymptotic form of the coefficients is conjectured to be 

~ A n n^ cf ] /n i/2 + (-4)V- 5fc - 9 £ df ] jn\ 

i>0 i>0 

where the amplitude estimates are given in Table f|. 

We have also determined the exact generating function for staircase polygons punc- 
tured by a single hole of perimeter 4, and also the generating function for staircase poly- 
gons punctured by a single hole of perimeter 6. We obtained these generating functions 
by generating the coefficients using the algorithm discussed, and then searching for an 
underlying differential equation. As a result we find the perimeter generating functions 
Pl(x) and Pq{x) for 1-punctured staircase polygons with a hole of perimeter 4 and 6 
respectively. 



and 



„„ , 2a; 4 - 16a; 3 + 20x 2 - 8x + 1 1 - Qx + 10x 2 - Ax 3 , , 

Pa(x) = ; ; , (11) 

41 ; 2(1 -Ax) 2VT^ 



„ s/ . 1 - 26x + 228x 2 - 906x 3 + 1709a; 4 - 1378x 5 + 322x 6 
PSW = 



2(1 -4a;) a 

32x 6 - 404a; 5 + 815x 4 - 586a; 3 + 182x 2 - 24a; + 1 
2(1 -Ax) 2 ' 



(12) 



Note that both these exact solutions display the confluent square root correction that we 
have found in our numerical investigations in the more general case. In the next subsection 
we analyse the analogous generating function for SAP. 



4.2 Self-avoiding polygons 



For unpunctured SAP, the perimeter generating function was recently extended |17j to 90 
step polygons, and the asymptotics clearly identified. The polygon generating function is 
defined to be 

P(°)( X ) = J2P2n% n ~ ^ (0) (^) + £ (0) 0*0(1 - /iV)i, (13) 
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where the functions and are believed to be regular in the vicinity of x c = 1/fJ?. 
We estimated fi = 2.63815853034(10). From this equation follows the asymptotic form of 
the coefficients, 

P?n ~ /i 2n n"t [ci + c 2 /n + c 3 /n 2 + c 4 /n 3 + ••■]. (14) 

We show in Table ^| the sequence of estimates of q, and in Table f| our estimates of the 
amplitudes, being the limits of the sequences {q}. 

The area generating function was first studied in ||, where the first 20 terms were 
given, and the asymptotic form estimated to be 

A<-°\q) = J> n <f ~ D(q) + E(q) log(l - nq), (15) 

n 

where k ~ 3.97087, and the logarithm in the above equation was understood to include 
the possibility of some power of a logarithm other than unity. (Though our analysis below 
implies that this is not the case.) In the present work we extend the series to 42 terms. 

4.2.1 Area generating function 

Using our greatly extended 42 term series, our analysis of the unpunctured disc area gen- 
erating function was carried out by standard methods. We used differential approximants 
|T3"| and found unbiased critical point and critical exponent estimates. The unbiased ex- 
ponent estimate had absolute value less than 10~ 5 , totally supporting our view that it is 
exactly zero. Assuming this, a biased estimate of the critical point is possible, and in this 
way we estimate « = 3.97094397(9). 

We then proceeded to seek the asymptotic form of the coefficients by writing 

a n = [q n }A(°\q)~ K «/nY,Ci/n m - 

i>0 

Our numerical results were well converged, demonstrating very convincingly that 
f(i) = i. This is the asymptotic form consistent with a pure logarithmic singularity, 
not raised to any power. Estimates of the amplitudes are given in Table [|. 

For 1-punctured discs, our analysis, based on the series known to q 45 , convincingly 
suggests the following asymptotic form: 

A^(q)^D^(q) + E^(q)/(l-Kq), 

and hence 

a£M? B ]^)~« n I>/n\ 

i>0 

Attempts to fit to alternative forms, corresponding to a confluent logarithm or a non- 
analytic correction-to-scaling term were unsuccessful, adding to our confidence that the 
above form is correct. Estimates of the amplitudes are given in Table |j. 

For twice-punctured discs, our analysis, also based on the series known to g 48 convinc- 
ingly suggests the following asymptotic form: 

A (2 \q)^ D ^(q) + E^\q)/(l-Kq) 2 , 
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and hence 

a® = [q n ]A^ 2 \q) ~ « n n^Q/n\ 

The conjectured asymptotic form for k punctured polygons, by area, is thus 

A {k \q) ~ DW(q) + E^(q)/(1 - nqf . 

4.2.2 Perimeter generating function 

The results for unpunctured polygons are fully discussed at the beginning of subsection 
[Ol As we found with staircase polygons, the generating function for punctured discs by 
perimeter was a more challenging numerical analysis problem than either its unpunctured 
counterpart, or its area counterpart. 

We found that the method of differential approximants was not particularly satis- 
factory. Given that we needed some 100 terms to successfully analyse the (presumably 
simpler) problem of punctured staircase polygons by perimeter, it is not surprising that for 
punctured SAP, for which we have 33 non-zero coefficients (corresponding to perimeters 
up to 84 steps), the method was not satisfactory. However, it did indicate the presence 
of a confluent singularity. As we found a confluent square-root singularity for staircase 
polygons, it is hardly surprising that a confluent singularity is detected for the punctured 
SAP generating function. In fact an exponent shift of around 1.7 was seen, compared 
to the expected value 1.5. We attribute this to the "short" series, coupled with the well 
known deleterious effect of confluent terms in such an analysis. Nevertheless, subsequent 
analysis of the asymptotic form of the coefficients, assuming an exponent shift of 1.5, 
together with a square-root confluent term, as found for punctured staircase polygons, 
gave satisfactory results. 

We denote the generating function for /c-punctured SAP, by perimeter, as 

p( fe )( x ) = J^pS^ ~ B (k) (x) + C {k \x)(l - f i 2 x) 1 - 5 - 1 - 5k , (16) 

n 

where the exponent is conjectured. For 1-punctured SAP, the vanishing of the expo- 
nent implies a logarithmic singularity. We fitted the coefficients to the asymptotic form 
appropriate to log(l — fi 2 x), so that the asymptotic form of the coefficients just involves 
decreasing integer powers of n. We then found that the estimates of the leading amplitude 
were monotonically increasing, which implies that the asymptotic form is wrong — too 
weak. Including a confluent square root singularity, as was found for punctured staircase 
polygons, stabilised the estimates. Accordingly, we conjecture that the asymptotic form 
is dominated by a logarithmic singularity, with a sub-dominant square root singularity, 
so that 

PS = [x n ]P {1 \x) ~ // 2 7n^>/nj (17) 

i>0 

Estimates of the amplitudes q are given in Table |4j. 

For twice punctured discs, a similar analysis suggested that the asymptotic form of 
the generating function is 

P( 2 )( X ) = J2P?nX n ~ B {2 \x) + C (2) (X)(1 - ^Vp 1 - 5 , (18) 
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again with evidence of a square root confluent term. As for 1-punctured SAP, we give 
estimates of the amplitudes q defined by 

in Table |. 

For thrice punctured discs, a similar analysis suggested that the asymptotic form of 
the generating function is 

P(3)( X ) = J2p ( 2nX n ~ B^ix) + C i3 \x)(l - ^ 2 XY\ (19) 
n 

again with evidence of a square root confluent term. As for 1-punctured SAP, we give 
estimates of the amplitudes q defined by 

i>0 

in Table |. 

The conjectured asymptotic form for k punctured polygons, by perimeter, is then 

p( k \ x ) = J2p { £x n ~ B {k \x) + C {k \x){l - fa) 1 - 5 - 1 ' 5 * (20) 

n 

where for k > we find strong evidence for a square root correction-to-scaling term. 



4.3 Polyominoes 

The problem of polyominoes has a long and interesting history, and has been well discussed 
in the popular scientific literature fL0| 



The enumeration of square lattice polyominoes to 24 steps [^3] was given in 1981 



extended to 25 steps in 1995 j3J and currently stands at 28 steps |f20[|. 

We have analysed the latest series by the method of differential approximants, and 
find the generating function behaves as 

V(y) = a n y n ~ G(y) + H(y) log(l - ry), (21) 

n 

where a n is the number of polyominoes of area n. In Q the estimate r = 4.06265(5) 
was given. The extra terms now available allow us to make the refined estimate r = 
4.062591(9). Analysis of the asymptotic form of the coefficients is totally consistent with 
a simple logarithm in the generating function. Thus 



[y n ]v(y)-r n J2^/ n 



i+1 

i>0 



The amplitude estimates are given in Table |4j. The leading amplitude is in complete 
agreement with, but 3 orders of magnitude more accurate than that given in ||TTJ, while 
the order of the leading term — 0(l/n) — was predicted by physical arguments back 
in 1981 [51]. In that work, the logarithmic singularity in the generating function of 
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strongly embedded site-animals was obtained. As discussed in the introduction, these are 
isomorphic to polyominoes. 

As well as extending the polyomino series, Oliveira e Silva has enumerated fc-punctured 
polyominoes for k < 6. Clearly, k = polyominoes are just SAP, and as we have 
seen in the previous section, these grow as K n where k = 3.9709.. < r. The arguments in 



25], |26| can be modified and applied to show that a fc-punctured polyomino has the same 
growth constant A as its unpunctured counterpart, and this is done in the Appendix. 
Thus the situation is that, for any finite number of punctures, 

but that 



\y n }V{y) = £ of) ~ r n . 



k>0 



We have analysed the series for £;-punctured polyominoes, k = 0,1, 2, and find the 
asymptotic form of the coefficients to be 



= {y n ]^ k \y) ^ K n n k -^^ ] In 1 



i>0 

,th 



corresponding to a generating function for /c-punctured polyominoes having a k order 
pole, viz: 

V^(y) = ai n ] y n ~ H^ k \y){l - Ky)-\ (22) 

n 

where k = is to be interpreted as a logarithm. Thus just as for punctured polygons, it 
is found that the exponent of fc-punctured polyominoes increases by 1 for each puncture. 
This is also proved in the Appendix. 

The coefficients q are given in Table 01. 



5 Conclusion 

We have investigated the effect of punctures on SAP and staircase polygons enumerated 
both by area and perimeter. 

In order to do this, we have developed a new algorithm, exponentially faster than 
direct counting, whereby we have radically extended a number of series. This extension 
was necessary in order to probe some rather subtle numerical behaviour. 

We found that, in every case, a finite number of punctures does not change the ex- 
ponential growth factor associated with the unpunctured counterpart of the punctured 
object being enumerated. 

This latter conclusion was also proved for polyominoes. The effect of punctures was 
also investigated numerically for polyominoes. 

Writing bn for the n th coefficient in the generating function for some k— punctured 
object, so that 

we found = tu^ in all cases. We found, further, that 6{k) = 9(0) + k if enumerating 
any of the objects we have considered by area. This can be proved, though subject, in 
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some cases, to the existence of the exponent in question. Subject to certain unproved 
assumptions we also showed that 8(k) = 9(0) + 3k/2 if enumerating by perimeter. 

We have, for the first time, obtained good numerical estimates of the sub-dominant 
terms for a range of problems, thus identifying both the nature of the generating function 
and any correction-to-scaling terms. 

We have also obtained an exact solution for the generating functions of staircase 
polygons, enumerated by perimeter, punctured by a single hole of perimeter 4 and of 
perimeter 6. These exact solutions provide additional support for our numerically based 
conjectures of the correction-to-scaling exponent in the general case. 

A more accurate estimate of the growth constant for SAP enumerated by area has 
been given, complementing our earlier work on the perimeter growth constants ||17|| . 



E-mail or WWW retrieval of series 

The series for the various generating functions studied in this paper can be obtained via 
e-mail by sending a request to I.Jensen@ms.unimelb.edu.au or via the world wide web on 
the URL [http: / /www. ms.unimelb.edu.au/ ^iwan/| by following the instructions. 
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Appendix. Growth constants and exponents for k- 
punctured polyominoes 

In this appendix we show that the growth constants are the same for all fc-punctured 
polyominoes. Our method of proof is based on that of van Rensburg and Whittington 
p5| , |26[| , making the necessary changes for the polyomino problem, and discussing in 
detail certain special cases. That is to say, if k denotes the connective constant for SAP 
enumerated by area, then this is the growth constant for fc-punctured polyominoes for any 
finite k. Further, if the usual asymptotic form for the number of fc-punctured polyominoes 
is assumed, sj? ~ CkrT^ h K n , and 0o exists, then <pk = 4>o — k. 



A.l Operations and mappings on punctured polyominoes 

(k) 

Let the set of all fc-punctured polyominoes with n cells be denoted by $„ , and the set 

(k) (k) 

of all polyominoes of n cells be $ n . Then, for each n, <3> n = \j k $„ . Throughout, let s n 
denote the cardinality of $n and let s n denote the cardinality of $ n . Hence s n = J2k s ^ ■ 
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Following [25|, we now define operations on punctured polyominoes, including con- 
catenation, drilling and surgery, which are needed in subsequent proofs. Concatenation 
allows us to change the size of polyominoes, while drilling and surgery are concerned with 
changing the number of holes in a polyomino. 

A. 1.1 Concatenation 

The concatenation mapping defined here is similar to that in ||14| . Consider the bounding 
rectangle R{P) of any polyomino P G $i . Define the top (bottom) edge of P to be the 
top (bottom) edge along the east (west) side of P(P). 

Now, the concatenation of two polyominoes P G $n and Q G &m is defined by 
joining P and Q while superimposing the top edge of P and the bottom edge of Q. The 
result is an (h + £;)-punctured polyomino with m + n cells, see figure [Al[ Hence we have 
a map 

T : $W x i — ► (Al) 

n m m+n \ J 

Lemma Al For all non-negative values of h and k, Sn Sm < s^+„ and in particular, 

„C0 > >) 

b n+l — bn 

Proof: For the mapping T defined above, every pair of polyominoes in the domain can 
be concatenated to form a larger polyomino in the codomain. Conversely, every such 
polyomino in the codomain can be uniquely broken up into the original ones. However, 
there are some polyominoes in the codomain which cannot be formed by concatenating 
two smaller polyominoes. An example is a 2 x 2 polyomino. Hence we get the first part 
of lemma as 

I n ml — l^ra+n r 

Putting k = and m = 1 and noting that |$^| = 1, we get the second part. ■ 
A. 1.2 Drilling 

Simply creating a hole inside an unpunctured polyomino by removing some interior cells 
does not allow us to drill certain classes of polyominoes. For instance, a polyomino 
composed of a single linear sequence of cells cannot be 'drilled' since removing any cell 
will either disconnect the polyomino or shorten the sequence. 

The following definition of drilling differs somewhat from that in [25L EB|, though the 
underlying idea is the same. 

First, we drill one hole. Let Po G be the polyomino that we are to puncture. 
Cover Po by a grid system, with each grid square of size 6x6 cells for any 6 > 5 (the 
minimum size will be justified later). Say 6 = 5. Pick any grid square G that covers at 
least one cell of Po and drill a hole there, as detailed below. 

Step 1 Remove all cells within the selected grid square G. If, after this step, we are 
left with a 1-punctured polyomino, we are done. If not, go to step 2. 

Step 2 Check the corners of G: If there are 2 disconnected components touching each 
other only at the corner, connect them by adding a cell at the appropriate corner of G. 
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Step 3 Put a 1-punctured polyomino with 8 cells at the centre of G. 

Step 4 Reconnect the disconnected components outside G to the 1-punctured poly- 
omino by adding linear sequences of cells (non- unique) . These steps are illustrated in 
figure 



Note that the minimum value of b is 5 because if we have anything less than that, 
we might create extra holes unexpectedly as the 1-punctured polyomino in Step 3 must 
touch the boundary of the grid square. An example in figure |A3| will illustrate this. In 
this example the original polyomino has no holes and the drilled polyomino has two holes, 
whereas with a grid of size 5 (or more), the number of holes in the drilled polyomino is 
only one. 

After this operation, the maximum number of cells that could be removed is b 2 = 25 
(finishing at step 1). The maximum number of cells that could be added is b 2 — 2 (this 
occurs when there was only 1 cell in the square before drilling and we end up with a 
1-punctured polyomino with b 2 — 1 cells after the operation). So, depending on each 
instance of the drilling operation, we obtain a resulting polyomino P\ G $j where j 
could be anything from n — b 2 to n + (b 2 — 2). The drilling operation thus defines a map 



-(6 2 -2) 



D : $1°) 



U 



(i) 



(A2) 



-6 2 



Theorem Al There exists a real constant C such that for all b > 5, 

< Cs^ +(62 _ 2) . (A3) 

Proof: Consider the intermediate (possibly disconnected) polyomino after Step 1, 
call it P l . There are many possible initial polyominoes which gives the same P l , i.e., all 
polyominoes with the same configuration outside G. They all map to "almost" the same 
resulting polyomino (the uncertainty implied in "almost" comes from the non-uniqueness 
of reconnecting in Step 4 which will be discussed later). Therefore, the mapping from 
domain to codomain is M'-to-one where M' is bounded by the number of ways that at 
most b 2 cells can be connected within G. Let the bound be M. 

On the other hand, Step 4 of the drilling process is not unique. There might be more 
than one way to connect those disconnected components to the punctured polyomino. 
But since G is finite, the number of ways of reconnection is bounded above. Let that 
bound be C' . 

Together with the drilling mapping D defined above, we can write, 

n+(fe 2 -2) 

s^KC'M £ sf. 

j=n—b 2 

By the increasing monotonicity of Sn over n (proven in lemma (|AT|)), 

S f < (2b 2 - l)C'Ms%,_ 2) . 
Set C = (2b 2 - l)C'M and the result follows. ■ 
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Now consider the drilling of h holes. Denote by |~c] the ceiling of c, the smallest 
integer greater than or equal to c. Now, we could choose h drilling locations from at least 
\n/b 2 ] grid squares. To see this, consider first a polyomino with n cells, where n < b 2 . 
The situation where we have the least number of drilling sites is when n cells fall exactly 
within 1 grid square. In this way, we only have 1 possible grid square where we could 
drill holes. Similarly, if we have b 2 n + w cells, where n > 1 and < w < b 2 , the minimal 
number of grid squares is when b 2 n cells fall into exactly n grid squares and the other w 
cells falls into another single grid. Then we have n + 1 possible grid squares to drill holes. 
Therefore we have at least \n/b 2 ] grid squares where we could drill holes. 

Letting H be a set of grid squares, let <£>«, denote the set of all polyominoes with n 
cells and \H\ holes where there is a hole in each grid square of H. Let s{f^ denote its 
cardinality. One property of this set is that, for any n\,n 2 G Z + such that n\ < 112, 



s [H] < s [H] 

Theorem A2 There exists a real constant K such that for all b > 5, 

( r "f < *4V* v»<r»/n. (A4) 

Proof: Place the polyomino in the grid system. Let A' be the set of all grid squares where 
we could drill holes. From previous results, we know there are at least \n/b 2 ~\ elements 
in A'. Truncate the set A' with only the first \n/b 2 ] elements and call this set of drilling 
sites A. 

Pick a subset H C A such that \H\ — h (< \n/b 2 ~\) and drill holes one by one in each 
grid square in H, leading to a series of mappings: 

M : $(°) ^ ^ $(2) ^ ... ^ 

where m/s are appropriate constants depending on each instance of operation and m« < 
n + (b 2 — 2)h, Vi. So, from this mapping M, we have 



Hence 



E 


s (0) 








s (0) 




b n 




M 


h ) 





S n+(b ,2 -2)h 



E[H] 
S n+(b 2 -2)h 



- - L * J n+(b 2 -2)h- 

The last line arises since {J HeA <$>f ] C = $^ and put K = C h . 
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A. 1.3 Surgery 



The surgery operation removes a linear sequence of cells inside the polyomino and con- 
catenates the sequence to the external boundary of the polyomino. Our objective is to 
join two holes thereby reducing the number of holes by one. 
We divide our domain into three classes: 

(1) the set of polyominoes which have at least one hole with size one cell; denote this 
set by with cardinality si, see figure |A"4](a); 

(2) the set of polyominoes which are not in (1), and have three holes touching each 
other at corners around a single cell; denote this set by $n with cardinality s„ , see 
figure 0b); 

(3) the set of polyominoes which are not in (1) nor (2) (the general case), denote this 
set by = U$i fe) ) with cardinality a#°. 

The general case: First, let's look at polyominoes in Consider a polyomino 

dn^ G ■ Let L be the set of all loop- free sequences^] of cells in a n h \ one end of which 
must touch the boundary of one hole and the other end must touch the boundary of 
another hole (or the exterior boundary). Define the length of the sequence to be the 
number of cells in the sequence. Let the set Z be the set of sequences in L that has 
minimum length. 

Now pick one sequence z G Z and cut it out. This step will join two holes together 
(or join one hole and the exterior together). Next, using the definition in section [A.l.l| , 
concatenate \ cz © D © z (D is a 3 x 3 polyomino block), and the resulting polyomino 
has one less hole and nine more cells than the original one. So we have a mapping 

S 3 : $00 ,_> $^-1) 

^ n n+9 

This mapping is at most ri-to-one because for any polyomino in the codomain, we could 
find at most n locations to connect z back. As a result, we find that, 

3? ) <™&? ) - (A6) 



Special cases: There are some subclasses of case (1) where we cannot apply the general 
surgery operation. One example is shown in figure ~K5 [27]. In this example, the minimum 



length of the loop-free sequence connecting two boundaries is two, but when we try to 
remove any such minimum sequence, we will disconnect the polyomino. 

To deal with case (1) polyominoes, we simply fill such a hole with a single cell. This 
satisfies our objective of reducing the number of holes by one. So we find a mapping 

Sl ; $00 | > ^ k - 1 \ 

KJ 1 n ^n+1 

The mapping Si is n' to 1 where n' is less than n. To see this, consider a polyomino in the 
codomain. To map back to the domain, we can choose any interior cell to remove and 
the number of choices is obviously less than the total number of cells in the polyomino, 
so n' < n. Therefore, 

s^<nst l) - (AT) 



1 A loop- free sequence is a sequence of cells where each successive pair of cells are joined and no cell 
appears more than once in the whole sequence 
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Finally, consider polyominoes in . The problem with these polyominoes is that 
when we try to remove the cell between the 3 holes, we will inevitably join 3 holes together 
instead of joining 2. One way to get around this is to deliberately choose another sequence 
(which is also a cell in this class) in the set Z. In particular, we choose to remove the 
cell z' G Z such that z' is the top cell (in lexicographic ordering) of Z which is not the 
problem cell. Similarly to the general case, we find a mapping 



S 2 ; $M 

* n n+9 



and hence 



Adding (TA"6|) , (|A7]), fl£|), we have 



and hence: 



Theorem A3 T/ie number of polyominoes of area n with h holes is bounded above by 3n 
times the number of polyominoes of area (n+9) with one less hole. That is, 

«? } < (A9) 
A. 2 Growth constants (by area) 

The following theorem proves the existence and equality of all the growth constants for 
/c-punctured polyominoes for all finite k. 

Theorem A4 There exists a constant (3$ such that for all h > 

lim n" 1 log(4 h) ) = log(A,). (A10) 



Proof: We use induction. First, (3q exists [p^| . Assume fit = lim^oo n 1 log(s n ) exists. 



From the results of the concatenation and surgery operations, we have 

,00 ,(1) < „(A+1) < o W 
a n-rn a m - s n — d/w, n+9- 

Choose some value m such that < Sm < oo, for example m = 8. Take the logarithm, 
divide by n and take the limit n — > oo. This gives 

log(/3 h ) < log(^ +1 ) < log(/3 A ) 

and hence /3ft = Iterating from j3 gives /3 = Ph for all finite h. ■ 
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A. 3 Relationships between critical exponents 

The following theorem establishes the relationships between critical exponents should they 
exist. 

Theorem A5 Assume for all h s [ n ] ~ C h n-+»ffi where C h is a h- dependent constant. 
Then 

<ph = 0o - h. (All) 
Proof: From the results of the drilling operation, we have 



Since 



h J h\ 



substituting the assumed asymptotic form Sn ~ C^n ^ h f3^ into (|A12|) , dividing by (3$, 



taking logarithms, dividing by log(n), letting n — ► oo and using the above result, we get 

h - O < -<f>h- (A13) 

Next, from the result of the surgery operation, we have 

4 h) < Zns% 1] < (3n)3(n + 9)^-1? <••■ 

< (3n)3(n + 9)...3(n + ^8%^ < [3(n + 9h)} h s^ +9h . 

Again by substituting the assumed asymptotic form, dividing by (3^\ taking logarithms, 
dividing by log(ra) and letting n — > oo, we get 

-<Ph<h- O . (A14) 

Hence 4>h — 0o — h. ■ 
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Table 1: The various 'input' states and the 'output' states (with corresponding weights) 
which arise as the boundary line is moved in order to include one more vertex of the 
lattice. 



Input 


Outputs 




'00' 


'00' 


:r 2 '12 


'Ol'/'IO' 


x'01' 


x'10' 


'02'/'20' 


x'02' 


x'20' 


'll'/'22' 


'00' 




'21' 


'00' 




'12' 


'accumulate' 





Table 2: A fit to the asymptotic form p 2 „ ~ /i 2 n 2 [c + ci/n + c 2 /n 2 + c 3 /n 3 + • • ■ ] for 
the number of SAP enumerated by perimeter. Estimates of the amplitudes Co, c±, c 2 , C3. 



n 




Co 




Cl 




c 2 




C3 


20 





09940085 


-0 


02745705 





02476376 





11822181 


21 





09940118 


-0 


02747548 





02511347 





11601107 


27 





09940177 


-0 


02751355 





02593211 





11011880 


28 





09940179 


-0 


02751510 





02597236 





10977030 


29 





09940180 


-0 


02751619 





02600168 





10950667 


30 





09940181 


-0 


02751694 





02602273 





10931043 


31 





09940182 


-0 


02751745 





02603734 





10916929 


32 





09940182 


-0 


02751777 





02604692 





10907354 


33 





09940182 


-0 


02751795 





02605254 





10901552 


34 





09940182 


-0 


02751802 





02605500 





10898929 


35 





09940182 


-0 


02751802 





02605494 





10898993 


36 





09940182 


-0 


02751796 





02605285 





10901358 


37 





09940182 


-0 


02751785 





02604913 





10905699 


38 





09940182 


-0 


02751771 





02604408 





10911757 


39 





09940182 


-0 


02751755 





02603796 





10919302 


40 





09940182 


-0 


02751736 





02603097 





10928158 


41 





09940182 


-0 


02751717 





02602327 





10938160 


42 





09940181 


-0 


02751696 





02601500 





10949174 


43 





09940181 


-0 


02751675 





02600629 





10961079 


44 





09940181 


-0 


02751653 





02599720 





10973796 


45 





09940181 


-0 


02751631 





02598785 





10987195 
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Table 3: A fit to the asymptotic form d$ ~ r] n n[co + Ci/nh + C2/n + C3/77J + c^/n 2 + • • • ] 
for the number of 1-punctured staircase polygons enumerated by area. Estimates of the 
amplitudes c , Ci, c 2 , c 3 , c 4 . 





Cn 






Ci 




Co 








C/i 




94 


6 83279x10" 


-3 


-1 


86263x10- 


-2 


-2 63249x10- 


-2 


2 95964x10- 


-2 


8 60593x10" 


-2 


95 


6.83258x10" 


-3 


-1 


86182x10" 


-2 


-2.64419x10" 


-2 


3.03465x10- 


-2 


8.42557x10" 


-2 


96 


6.83256x10" 


-3 


-1 


86173x10" 


-2 


-2.64559x10" 


-2 


3.04363x10- 


-2 


8.40388x10" 


-2 


97 


6.83239x10" 


-3 


-1 


86109x10" 


-2 


-2.65488x10" 


-2 


3.10389x10- 


-2 


8.25743x10" 


-2 


98 


6.83235x10" 


-3 


-1 


86092x10" 


-2 


-2.65738x10" 


-2 


3.12014x10- 


-2 


8.21773x10" 


-2 


99 


6.83222x10" 


-3 


-1 


86040x10" 


-2 


-2.66501x10" 


-2 


3.17010x10- 


-2 


8.09504x10" 


-2 


100 


6.83217x10" 


-3 


-1 


86019x10" 


-2 


-2.66815x10" 


-2 


3.19076x10- 


-2 


8.04405x10" 


-2 


101 


6.83206x10" 


-3 


-1 


85976x10" 


-2 


-2.67457x10- 


-2 


3.23325x10- 


-2 


7.93864x10" 


-2 


102 


6.83200x10" 


-3 


-1 


85953x10" 


-2 


-2.67801x10- 


-2 


3.25612x10" 


-2 


7.88161x10" 


-2 


103 


6.83191x10" 


-3 


-1 


85916x10" 


-2 


-2.68355x10- 


-2 


3.29313x10- 


-2 


7.78885x10" 


-2 


104 


6.83185x10" 


-3 


-1 


85893x10" 


-2 


-2.68708x10- 


-2 


3.31690x10- 


-2 


7.72899x10" 


-2 


105 


6.83177x10" 


-3 


-1 


85860x10" 


-2 


-2.69197x10- 


-2 


3.34985x10- 


-2 


7.64560x10" 


-2 


106 


6.83171x10" 


-3 


-1 


85837x10" 


-2 


-2.69553x10- 


-2 


3.37400x10" 


-2 


7.58417x10" 


-2 


107 


6.83164x10" 


-3 


-1 


85809x10" 


-2 


-2.69988x10" 


-2 


3.40364x10" 


-2 


7.50842x10" 


-2 


108 


6.83159x10" 


-3 


-1 


85786x10" 


-2 


-2.70333x10" 


-2 


3.42727x10- 


-2 


7.44775x10" 


-2 


109 


6.83153x10" 


-3 


-1 


85761x10" 


-2 


-2.70731x10- 


-2 


3.45466x10" 


-2 


7.37709x10" 


-2 


110 


6.83148x10" 


-3 


-1 


85739x10" 


-2 


-2.71063x10" 


-2 


3.47758x10" 


-2 


7.31769x10" 


-2 


111 


6.83142x10" 


-3 


-1 


85716x10" 


-2 


-2.71427x10- 


-2 


3.50291x10- 


-2 


7.25171x10" 


-2 


112 


6.83137x10" 


-3 


-1 


85696x10" 


-2 


-2.71747x10- 


-2 


3.52522x10- 


-2 


7.19335x10" 


-2 


113 


6.83132x10" 


-3 


-1 


85675x10" 


-2 


-2.72077x10- 


-2 


3.54832x10- 


-2 


7.13266x10" 


-2 
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Table 4: Amplitude estimates appearing in the asymptotic form of the coefficients of the 
various models considered. The order of the term associated with the amplitude q varies 
from model to model, and is given in the text for each model. The various connective 
constants are: rj = 2.30913859330, re = 3.97094397, n = 2.63815853034, r = 4.062591. All 
numbers quoted are expected to have errors only in the last quoted digit, k is the number 
of punctures. 

The n th coefficient is given by Prefactor x ^2 i>0 Ci/n^ . See the text for the problem 
dependent value of f(i). 



Model, 


k 


Prefactor 


c 


Cl 


c 2 


c 3 


Parameter. 














Staircase 







0.12881579 


0(1.5996 n ) 






polygons 


1 


r) n n 


0.006831 


-0.0185 


-0.028 


0.04 


by area 


2 


r] n n 2 


7.87 x 10- 5 


-0.00043 


0.0010 


-0.007 




3 


rf n 


6.04 x 10" 7 


-5.03 x 10~ 6 


0.000031 


-0.00022 


Staircase 





4 n n-i/v^r 


1/4 


3/32 


25/512 


64/4096 


polygons 


1 


4™ 


0.0147 


-0.19 


1.4 




by 


2 


4« n f 


8.0 x 10~ 4 


-2 x 10~ 2 


0.2 




perimeter 


3 


4 n n 3 


3.0 x 10~ 6 


-1 x 10~ 3 


0.02 




SAP 







0.408105 


-0.5467 


0.626 


-3 


by 


1 




0.000975 


-0.0097 


-0.04 




area 


2 




0.00000118 


-0.000019 


-0.0001 






3 


K n n 2 


1.0 x 10~ 9 


-2 x 10~ 8 


-3 x 10~ 7 




SAP 





2n — — 

fx zn n 2 


0.0994018 


-0.02751 


0.0255 
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Figure 1: An example of a punctured self-avoiding polygon (left panel) and a punctured 
staircase polygon (right panel). The thick lines show the perimeter of the enclosing 
polygon while the medium lines show the perimeter of the holes. 



Figure 2: An example of a punctured polyomino that is not a punctured polygon 



Figure 3: A snapshot of the intersection (dashed line) during the transfer matrix cal- 
culation on the square lattice. Polygons are enumerated by successive moves of the kink 
in the intersection, as exemplified by the position given by the dotted line, so that one 
vertex at a time is added to the rectangle. To the left of the intersection we have drawn 
an example of a partially completed polygon. 
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Figure 4: Some of the local configurations which occur as the kink in the intersection is 
moved one step. 




Figure 5: Illustration of how a pair of loops can be placed (left), connected to produce 
a valid SAP (middle), and connections leading to forbidden graphs (right). 
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Figure 6: Illustration of how connecting loop ends lead to various cases of punctured self- 
avoiding polygons. In the upper (lower) panels we show how the configuration {11112222} 
({1111212222}) can be completed given that the first pair of '12' edges are connected to 
form a puncture (in the lower panel we further connect the two 1-edges on either side of 
the puncture). In each panel we also show the intersection line with the numbers giving 
the labelling of the loop ends. 
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Figure Al: Concatenation operation 



(a) the original polyomino 



(b) after adding cells in step 2 



(c) after adding a 1-punctured polyomino (d) the final drilled polyomino after reconnections 



Figure A2: Drilling operation: (a) The original polyomino. After step 1, all the cells 
within the grid square are removed. In step 2, we reconnect all the disconnected com- 
ponents around the corner (shown in (b)). Then, in step 3, a 1-punctured polyomino is 
placed at the centre of the grid square (shown in (c)). Finally, all the disconnected pieces 
outside the grid square are reconnected to the punctured polyomino (shown in (d)). 



35 



(a) the original polyomino (b) the drilled polyomino 



Figure A3: An illustration of a grid square that is too small. 



(a) a polyomino in (b) a polyomino in $ 



(h) 



Figure A4: An illustration of two special cases for the surgery operation. In (b), P 
denotes the problem cell whose removal will join 3 holes together. In this case we choose 
to remove z' 



Figure A5: A polyomino in $^ where we cannot apply the general surgery operation 
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